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Abstract 

In  many  surveillance  and  tracking  systems  multi-sensor  configurations  are  used  to  provide 
a  greater  depth  of  information  on  the  target  and  likewise  to  increase  the  robustness  of  the 
sensor  network  to  survive  individual  sensor  failure. 

This  lecture  investigates  several  aspects  of  sensor  networks  and  illustrates  various  properties 
of  distributed  networks  composed  of  identical  sensors. 


1  Introduction 

Sensor  networks  are  very  attractive  for  reconnaissance  and  surveillance  applications  as  they 
provide  more  information  from  the  same  target,  which  helps  to  identify,  classify,  and  track 
it.  Furthermore,  sensor  networks  possess  a  higher  probability  to  survive  if  individual  sensors 
fail.  Looking  at  a  sensor  network  and  its  operation  these  can  be  divided  into  two  main 
categories: 

•  mono-type  and 
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•  multi-type  sensor  networks. 

Networks  composed  of  different  sensor  types  allow  to  determine  various  characteristics  of 
the  same  target  and  therefore  improve  the  classification  and  identification  task,  but  not  nec¬ 
essarily  the  target  localization  accuracy. 

Depending  on  the  arrangement  of  the  sensors  in  the  network  a  further  classification  can  be 
introduced  as 

•  colocated  network  and 

•  distributed  sensor  network. 

Any  combination  of  these  first  categories  is  possible,  but  generally  different  sensors  are 
located  nearby  to  observe  the  same  area.  An  example  of  a  colocated  mono-type  sensor  net¬ 
work  is  the  virtual  linear  array  antenna,  composed  of  several  transmit  and  receive  antennas 
arranged  along  a  straight  line  [1],  For  such  a  configuration  the  signal  processing  will  be 
performed  by  a  central  processing  unit,  which  ensures  that  the  network  operates  coherently. 
The  common  architecture  of  distributed  multi-sensor  networks  is  to  install  a  preprocessing 
at  each  receiver  node  and  transfer  the  results  over  a  high-speed  communication  link  to  a  fur¬ 
ther  processing  stage  for  data  fusion.  Depending  on  the  signal  bandwidth  and  complexity  of 
the  distributed  network  also  a  central  processing  solution  may  be  feasible,  depending  on  the 
communication  link  between  the  nodes  and  the  central  unit  [2] . 

The  advantage  of  a  distributed  mono-type  sensor  network  is  the  possibility  of  bi-/multi static 
operation.  That  implies,  a  receiver  node  is  able  to  receive  echoes  from  the  target  illuminated 
from  different  transmit  nodes.  Furthermore,  spatially  distributed  mono-type  sensor  networks 
offer  more  degrees  of  freedom  than  a  system  with  one  sensor.  These  multiple-input  multiple- 
output  (MIMO)  networks  support  flexible  time-energy  management  modes,  exhibit  higher 
angular  resolution,  show  an  improved  parameter  estimation,  and  facilitate  high  resolution 
target  localization  accuracy  with  an  improved  detectability  by  exploiting  target  echo  (RCS) 
diversity  [3].  An  optimal  MIMO  network  comprises  multiple  transmitters  and  multiple  re¬ 
ceivers  to  pick-up  the  reflected  signal.  Radar  systems  with  their  all-day,  all  weather  ability 
are  promising  sensors  for  MIMO  networks  for  surveillance  applications  in  defence  and  civil¬ 
ian  categories.  On  the  contrary,  if  a  networks  is  composed  of  different  sensor  types  each 
node  can  only  operate  in  a  monostatic  configuration. 

A  MIMO  network  can  be  split  up  into  two  principal  categories  of  operation  based  on  trans¬ 
mitter  usage  [10]: 

•  spatial  multiplexing  and 

•  space  time  transmit  diversity 
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In  the  spatial  multiplexing  case,  each  transmit  node  in  the  network  has  an  unique  orthogonal 
coded  signal,  which  is  transmitted  in  the  same  frequency  channel.  If  these  electromagnetic 
waves  arrive  at  the  receive  node,  the  receiver  can  separate  these  signals  by  matched  filters, 
creating  as  many  parallel  channels  as  transmit  nodes  exists.  Spatial  multiplexing  is  a  very 
powerful  technique  for  increasing  MIMO  network  capacity  and  can  be  used  with  or  without 
transmit  channel  knowledge. 

In  contrast  to  the  spatial  multiplexing  case  all  transmit  nodes  share  the  same  waveform  and 
frequency  channel  in  the  space  time  transmit  diversity  case.  Due  to  this,  always  only  one 
transmitter  is  active.  After  the  transmitted  and  reflected  signal  has  arrived  at  all  receive 
nodes  the  next  transmitter  will  emit  a  signal.  As  the  time  to  of  each  firing  is  known  the 
receivers  can  assign  the  received  signal  to  the  corresponding  transmitter  [3]. 

The  following  sections  of  this  lecture  will  focus  on  MIMO  radar  networks,  representative 
for  all  other  mono-type  sensor  networks,  and  will  discuss  some  essential  design  aspects  to 
obtain  the  optimal  performance  for  a  given  configuration. 


2  Synchronisation 

In  a  sensor  network  synchronisation  between  all  nodes  is  essential  for  the  performance  of  the 
system  and  can  be  parted  in  three  areas:  timing,  frequency,  and  phase.  Timing  synchronisa¬ 
tion  is  needed  between  transmit  and  receive  nodes  for  range  measurement  and  typically  an 
accuracy  in  the  order  of  a  fraction  of  the  transmitted  compressed  pulse  is  required  r  =  1/ B, 
where  B  is  the  signal  bandwidth  [4]. 

If  the  network  nodes  are  not  too  far  separated  from  each  other  synchronising  the  stable  local 
oscillators  can  be  achieved  by  connecting  them  together  by  cable,  fibre  or  direct  commu¬ 
nication  link.  If  there  is  no  direct  line  of  sight  available,  time  synchronisation  can  also  be 
achieved  via  a  scattering  signal  (clutter)  or  even  a  scattering  of  the  troposphere.  A  prereq¬ 
uisite  is  that  the  scattering  volume  is  detected  by  transmit  and  receive  antenna  lobes.  This 
method  is  not  suitable  to  stabilize  two  local  oscillators  via  a  phase-locked  loop  as  one  can 
aspect  a  large  variation  (jitter)  in  the  time  base  if  the  scattering  body  moves,  for  instance  by 
gusts  of  wind.  Commercially  available  are  a  range  of  various  qualities  of  stable  oscillators 
ranging  from  a  simple  quartz  oscillators,  temperature  controlled  quartz  with  single  or  double 
ovens  for  stable  operating  temperature,  to  more  expensive  atomic  clocks  such  as  rubidium  or 
caesium  oscillators. 

Due  to  the  inherent  aging  and  instability,  the  local  references  must  be  continuously  re¬ 
synchronised  on  a  time  interval,  which  depends  on  the  required  stability  and  coherence  in 
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the  network.  The  stability  directly  influences  the  coherent  integration  time  in  the  sensor 
network.  Another  important  aspect  of  reference  oscillators  in  a  distributed  MIMO  radar  net¬ 
work  is  their  phase  noise  characteristic.  While  in  a  single  sensor  system  and  its  monostatic 
operation  the  received  echo  has  some  correlation  with  the  transmitted  signal,  even  at  great 
distances,  such  a  correlation  does  not  exist  in  a  MIMO  network  with  several  independent 
reference  clocks. 

Phase  noise  is  a  representation  of  random  fluctuations  in  the  phase  of  a  waveform  and  is  spec¬ 
ified  in  frequency  domain  for  a  given  frequency  offset  from  the  carrier  e.g.  £  =  110  dBc  at 
10  Hz  and/or  in  time  domain  by  the  two-point  variance  (Allan  variance)  for  a  set  of  time  in¬ 
tervals.  Particularly  for  detecting  slow  moving  targets  which  introduce  a  small  Doppler  shift 
a  low  phase  noise  close  to  the  carrier  is  very  important  for  the  Doppler  processing.  Attention 
has  to  be  paid  when  the  system  is  mounted  on  a  vibrating  platform,  such  as  an  airplane,  as 
oscillators  are  sensitive  against  accelerations  and  phase  noise  will  degrade  [5].  A  comparison 
between  some  standard  oscillators  and  their  single  sideband  noise  (£  in  [dBc/Hz])  is  shown 
in  table  1 . 


Table  1:  Comparison  of  various  stable  oscillators  at  10  MHz 


TCXO 

ocxo 

OCXO(BVA) 

MCXO 

Rubidium 

Cesium 

GPS 

temperature 

stability 

MO-6 

2-10-8 

MO-10 

2-10-8 

2-10-8 

2-10-8 

drift  per  day 

1-10-8 

MO"10 

5T0-12 

5T0-11 

5T0-13 

3T0“14 

3T0-14 

short-time 
stability  ay{r) 

1-10“ 9 

1-10”12 

5T0-13 

l-lO"10 

5T0-12 

Is 

MO"10 

M0-11 

1-10-12 

2-10-8 

M0-11 

6T0-11 

— 

100  s 

3T0 1° 

l-lO"11 

l-lO-12 

2-10-8 

l-lO"13 

6T0“12 

— 

1  day 

1-10”12 

3T0“12 

3T0~12 

5TO-10 

l-lO”13 

6T0~12 

— 

£  @  1Hz 

-50 

-100 

-122 

-115 

-80 

-85 

— 

£  @  10Hz 

-80 

-130 

-137 

-135 

-98 

-125 

— 

£  @  100Hz 

-110 

-140 

-145 

-145 

-137 

-135 

— 

£  @  1kHz 

-120 

-145 

-156 

-150 

-150 

-140 

— 

If  no  direct  synchronisation  is  possible,  as  the  baseline  is  too  large  for  any  cable  or  fibre, 
or  no  line-of-sight  exists,  the  local  oscillators  of  each  network  node  can  be  indirectly  syn¬ 
chronised  by  use  of  a  Global  Navigation  Satellite  System  (GNSS)  which  provides  a  highly 
stable  time  reference.  This  time  reference,  the  1  Pulse-Per-Second  (PPS)  signal,  can  be  used 
to  discipline  the  local  oscillator.  The  uncertainty  in  time  of  a  1PPS  is  in  the  order  of  100  ns 
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or  less.  If  two  GNSS  receivers  are  receiving  navigation  signals  from  the  same  satellites, 
which  is  possible  up  to  a  distance  of  8000  km  between  them,  an  accuracy  of  5. ..20  ns  can 
be  achieved,  even  if  the  selective  available  (SA)  signal  is  activated.  Only  the  location  of  the 
GNSS  receiver  must  known  precisely  (<  1  m).  This  can  be  easily  obtained  by  an  averaging 
of  the  calculated  position  from  the  GNSS  data  over  several  days.  Figure  1  shows  the  phase 
difference  between  two  temperature  controlled  crystal  oscillators  which  were  disciplined  by 
different  GPS  receivers  and  their  PPS  signal.  The  standard  deviation  of  the  1PPS  signal  was 
50  ns  for  each  GPS  receiver.  Each  GPS  receiver  has  been  set  into  the  mobile  mode,  which 
means  that  the  receiver  has  to  solve  the  full  set  of  equations  to  determine  the  position  and 
time  for  each  GPS  epoch.  Therefore,  a  lower  variance  can  be  expected  if  both  GPS  units 
works  in  the  stationary  mode ,  where  the  location  of  the  GPS  receiver  is  known  and  only  the 
time  has  to  be  determined  from  the  GPS  signals  [6]. 


Figure  1:  Measured  phase  difference  between  two  GPS  disciplined  TCXO  at  10  MHz 


In  this  configuration  the  local  oscillator  determines  the  short  time  stability,  while  the  long 
time  stability  is  determined  by  the  stability  of  the  GPS  system.  As  long  as  the  receiver  is 
locked  to  the  GPS,  the  long  time  stability  is  ensured.  In  the  case  of  GPS  jamming  or  receiver 
unlock,  the  stability  of  the  local  oscillator  unit  will  determine  the  drift  due  to  aging  and 
temperature  changes  etc. 

2.1  Phase  synchronisation 

For  Doppler  or  MTI  processing  phase  coherence  between  transmit  and  receive  nodes  has  to 
be  established,  which  allows  to  reject  clutter  or  chaff.  The  phase  coherence  can  be  obtained 
in  the  same  ways  as  time  coherence.  Using  indirect  phase  synchronisation,  which  is  the  best 
solution,  involves  high-precision  oscillators  at  the  network  nodes  that  are  re-synchronised 
via  the  1PPS  provided  by  the  GNSS  receiver.  Over  the  whole  coherent  integration  time  T/,, 
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the  phase  stability  has  to  be  guaranteed  and  equal  to  A (j>  /  2Tcfrk  [4].  For  a  ground  based 
bistatic  radar  with  a  center  frequency  of  3  GHz,  a  maximum  phase  deviation  of  A©  =  4°  = 
0,  07  mrad  and  a  coherent  integration  time  of  A rk  =  1  ms  the  required  oscillators  stability  is 
3.7  •  10~12.  This  requirement  can  be  fulfilled  with  a  temperature-controlled  crystal  oscillator, 
as  can  be  taken  from  table  1. 

Likewise  the  phase  synchronisation  can  be  achieved  by  a  direct  signal,  if  a  line-of-sight 
exists.  If  the  synchronisation  takes  place  along  each  transmit  pulse  a  phase  stability  can 
be  reached  of  A©  /  2nfATrt,  with  A Trt  denoting  the  travelling  time  difference  between 
transmitter-target-receiver  and  the  direct  signal  transmitter-receiver.  A©  is  the  allowed  phase 
difference  in  rad.  If  an  accuracy  of  A©  =  4°  =  0,  07  rad  is  requested  at  a  center  frequency 
of  /  =  10  GHz  and  the  time  difference  A Trt  =  1  ms  (A rrt  =  300  km)  the  oscillator  must 
have  a  stability  of  10“9,  which  can  easily  be  achieved  by  a  simple  quartz  oscillator. 

3  Ambiguity  function 

An  important  tool  to  evaluate  radar  or  sonar  signal  processing  performance  in  terms  of  range 
and  Doppler  resolution  as  well  as  clutter  rejection  is  the  ambiguity  function.  The  concept  of 
ambiguity  function  was  first  introduced  by  Woodward  [7].  It  is  a  two  dimensional  function 
of  time  delay  and  Doppler  frequency  x(r,  fD))  showing  the  absolute  envelope  of  the  output 
of  the  receiver  matched  filter  when  the  input  to  the  filter  is  a  Doppler  shifted  version  of  the 
original  transmitted  signal.  The  ambiguity  function  is  determined  only  by  the  properties  of 
the  received  pulse  and  the  matched  filter,  which  represents  the  transmitted  pulse,  and  not  any 
specific  target  scenario.  Many  definitions  of  the  ambiguity  function  exist.  Several  focus  on 
narrowband  signals  and  other  are  applicable  to  describe  the  propagation  delay  and  Doppler 
relationship  of  wideband  signals.  For  a  complex  baseband  signal  s(t),  which  fulfils  the 
narrowband  condition  2  vBT /c  ■ Cl,  where  v  is  the  target  velocity,  B  is  the  signal  bandwidth, 
T  is  the  pulse  duration,  and  c  is  the  speed  of  propagation,  the  narrowband  ambiguity  function 
is  given  by  [7]: 


(1) 


where  *  denotes  the  complex  conjugate  and  fo  the  Doppler  shift  in  frequency.  An  important 
assumption  for  the  target  is  that  its  scattering  properties  do  not  change  over  the  pulse  duration 
and  with  the  look  angles  and  that  is  only  slowly  manoeuvring. 

The  monostatic  ambiguity  function  was  developed  for  a  single  colocated  transmit/receive 
pair  and  is  fairly  well  developed  and  understood  [8].  It  has  been  shown  that  the  ambiguity 
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function  arises  from  the  detection  and  parameter  estimation  problems  joined  with  a  slowly 
fluctuating  point  target  being  observed  in  additive  white  Gaussian  noise. 


For  bistatic  geometry  the  simple  relationship  between  time  delay  and  range  and  between  tar¬ 
get  velocity  and  Doppler  shift  is  no  longer  valid  and  therefore  Tsao  et  al.  [9]  have  developed 
the  bistatic  ambiguity  function.  It  is  shown  that  the  bistatic  geometry  has  a  strong  influence 
on  the  ambiguity  function  x(r,  fD),  which  can  be  written  as: 


X(RrHi  RRai  Vh,  Va,  Or,  L ) 


2 


Ra(RRA,  Or,  L)) 


S*(t  —  TH(Rrh,0r,  L)) 


■  exp{-j(uDH(RRH,  VH,  Or,  L) 


—uda(Rra,  Va,  Or,  L))  t}  dt 


2 


(2) 


in  which  Rr  and  Rj  are  the  ranges  from  the  target  to  the  receiver  respectively  to  the  trans¬ 
mitter,  V  is  the  target  radial  velocity,  0R  is  the  angle  of  the  target  measured  from  the  receiver, 
L  is  the  bistatic  baseline,  r  is  the  transmitter-target-receiver  delay  time,  and  the  subscription 
H  and  A  denote  the  hypothesized  and  actual  values. 


The  important  difference  between  the  monostatic  (1)  and  bistatic  ambiguity  function  (2)  is 
that  the  geometrical  layout  of  the  transmitter,  receiver  and  target  are  now  taken  into  account. 
This  can  have  a  significant  effect  on  the  form  of  the  ambiguity  function  and  the  resulting 
range  and  Doppler  resolutions.  For  the  following  examples  a  pulse  is  considered  with  a 
complex  envelope  described  by  a  simple  Gaussian  shape: 

s(i)  =  ^(-^2)  -oo<t<oo  ,  (3) 

with  T  the  pulse  duration.  The  monostatic  ambiguity  function  for  such  a  pulse  is  shown 
in  Fig.  2.  The  following  parameters  where  assumed  for  the  simulations  of  the  ambiguity 
function:  pulse  duration  T  =  120  /is,  time  delay  r  =  400  /is,  and  center  frequency  fc  = 
50  MHz.  All  figures  show  the  considered  geometry  (left  subplot),  the  normalized  ambiguity 
function  as  a  function  of  velocity  and  target  range  from  the  receiver  (center  plot),  and  an 
equal-height  contour  plot  (right  subplot).  As  the  linear  relationship  between  delay  and  range, 
and  Doppler  shift  and  radial  velocity  is  no  more  valid  in  the  bistatic  case  the  ambiguity 
function  is  plotted  with  respect  to  range  and  velocity.  For  the  monostatic  case  these  plots  are 
identical.  However,  when  targets  are  near  or  even  crossing  the  bistatic  baseline  the  bistatic 
ambiguity  function  changes  dramatically,  as  shown  in  Fig.  3  and  Fig.  4.  Especially  if  the 
target  crosses  the  baseline  perpendicularly  no  resolution  in  range  or  Doppler  can  be  obtained, 
hence,  the  bistatic  ambiguity  function  in  this  area  is  zero,  regardless  of  the  position  of  the 
target  on  the  baseline  and  the  velocity. 
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Figure  2:  Monostatic  ambiguity  function  plot  of  Gaussian  pulse:  left  figure  shows  the  geom¬ 
etry,  center  figure  a  three  dimensional  plot,  and  on  the  right  side  a  contour  plot  is  shown. 


Figure  3:  Bistatic  ambiguity  function  plot  of  Gaussian  pulse:  left  figure  shows  the  geometry, 
center  figure  a  3  dimensional  plot  in  the  range- velocity  plane,  and  on  the  right  side  a  contour 
plot  is  shown. 


For  a  MIMO  radar  network,  which  can  be  seen  as  a  composition  of  several  bistatic  transmit 
receive  pairs,  the  ambiguity  function  is  formulated  based  on  the  bistatic  ambiguity  function. 
It  is  assumed  that  the  network  is  composed  of  M  transmitters  and  N  receivers.  In  that  case 
the  network  shows  MN  bistatic  pairs.  To  simplify  the  derivation  of  the  multistatic  ambiguity 
function  the  same  assumptions  are  made  as  for  the  bistatic  ambiguity  function.  Furthermore, 
it  is  assumed  that  the  network  is  coherent.  This  implies  that  the  echoes  arriving  at  different 
time  instances  can  be  processed  jointly.  Similar  to  the  bistatic  radar  ambiguity  analysis,  the 
multistatic  radar  ambiguity  function  is  developed  by  the  following  three  steps  [10]: 

(i)  For  each  transmitter-receiver-pair  the  bistatic  ambiguity  (2)  function  is  calculated 

(ii)  Calculating  a  weighting  factor  according  to  received  signal  intensity 


,  _  Pi  Gi  Gj  A2 erg 

(47 

P 

1  11 

I  •  •  =  - - - 

II  Max(Pij) 


(4) 

(5) 
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Figure  4:  Bistatic  ambiguity  function  plot  of  Gaussian  pulse:  left  figure  shows  the  geometry, 
center  figure  a  3  dimensional  plot  in  the  range- velocity  plane,  and  on  the  right  side  a  contour 
plot  is  shown. 


(iii)  To  formulate  multistatic  radar  ambiguity  function  using  the  results  from  previous  cal¬ 
culations: 


Xmulti 


1 


M  N 


M 2  N 2 


Wi^ 


i= 1  j= 1 


(6) 


A  multistatic  radar  model  has  been  developed  for  a  comprehensive  understanding  of  multi¬ 
static  radar  ambiguity  performance.  It  is  assumed  that  the  fixed  transmitters  and  receivers 
are  located  in  one  plane  and  the  target  is  moving  in  another  plane,  which  is  parallel  to  the 
transmitter-receiver  plane.  This  tool  is  very  useful  to  evaluate  the  performance  of  a  given 
multistatic  network. 


4  Target  localisation 

Due  to  their  improved  target  parameter  estimation  capability  distributed  MIMO  radar  net¬ 
works  are  very  attractive  for  system  designers.  In  particular  their  improved  angular  resolution 
capability  and  ability  to  separate  multiple  targets  [12],  [13],  to  improve  parameter  identifica¬ 
tion  [14],  and  to  improve  radar  performance  by  exploiting  radar  cross  section  (RCS)  diversity 
is  very  attractive  for  surveillance  applications  [15].  Likewise  MIMO  radar  networks  can  han¬ 
dle  slow  moving  targets  by  exploiting  Doppler  estimates  from  multiple  directions  [16]  and 
feature  a  highly  accurate  estimation  of  target  position  [17],  [18]. 

The  estimation  of  target  localisation  in  a  MIMO  radar  network  can  be  performed  in  a  coher¬ 
ent  or  non-coherent  way  based  on  best  linear  unbiased  estimator  (BLUE). 
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The  distributed  MIMO  radar  network  consists  of  M  transmit  and  N  receive  nodes,  which  are 
located  in  the  two-dimensional  plane  (x,y).  The  transmitters  are  located  at  7).  =  (xtk,  ytk ), 
with  k  =  1, M,  and  the  coordinates  of  the  receivers  are  Ri  =  (xri,  yri ).  Each  transmitter 
emits  an  orthogonal  waveform,  with  a  lowpass  equivalent  sfc(f),  k  =  1,  All  wave¬ 

forms  are  narrowband  signals  fulfilling  the  assumption  of  (B  /  fc)2  -C  1.  The  target  at  the 
position  X  =  (x,  y)  has  a  complex  radar  cross  section  (,  which  does  not  change  over  the 
aspect  angle,  and  introduces  a  time  delay  r/fc  for  each  received  signal.  Taking  into  account 
these  assumptions  the  equivalent  lowpass  received  signal  at  receiver  l  can  be  described  by 
[20]: 

M 

ri(t)  =  y;  (sk(t  -  Tik)e~j2nfcTik  +  m(t)  (7) 

k= 1 


with  ri;  (f)  the  complex  Gaussian  white  noise.  The  target  location  can  be  computed  by  the 
time  difference  of  arrival  (TDOA)  of  the  received  signal  to  more  than  two  receivers.  In 
Eq.  (7)  the  time  delays  r/fc  ,  which  are  determined  by  the  location  of  the  target,  the  transmit¬ 
ters,  and  receivers,  are  given  by  the  following  relation: 


Tlk  ~  ~ 

C  L 


y/ (xtk  -  x)2  +  (ytk  -  y)2  +  yj {xri  -  x)2  +  ( yri  -  y)2 


(8) 


The  target  localisation  processing  uses  the  well-established  multilateration  principle  and  can 
be  performed  coherently  or  non-coherently.  In  the  non-coherent  case  each  receiver  node  in 
the  MIMO  network  has  only  timing  information  on  all  emitted  waveforms,  but  no  knowledge 
of  the  phase.  Hence,  as  no  phase  information  is  available  for  non-coherent  localisation  the 
received  signal  in  Eq.  (7)  simplifies  to: 


M 

ri(t )  =  ^2  aiksk(t  -  Tik)  +  ni(t)  (9) 

k= l 


where  au-  are  unknown  complex  amplitudes  effected  by  the  targets  radar  cross  section  and 
the  phase  offset  between  transmitter  and  receiver. 

The  goal  of  the  multi  radar  network  is  to  estimate  the  target  position  X  =  (x,  y).  This  can 
be  achieved  by  formulating  the  maximum  likelihood  estimate  (MLE)  using  the  equation  sets 
given  by  Eq.  (7).  Another  way  to  determine  the  target  location  is  to  estimate  first  the  time 
delays  r/^  from  the  equations  formed  by  Eq.  (8)  and  to  compute  in  a  second  step  the  position 
by  employing  multilateration  principles. 

Another  approach  proposed  by  Godrich  et  al.  [11]  uses  a  linear  model  in  estimating  the  time 
delays  rik  given  by  the  output  of  the  matched  filter  to  the  transmitted  signals  sk(t).  For  a 
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fully  coherently  operating  radar  network,  each  node  has  global  time  and  phase  information, 
and  the  estimation  is  given  by: 

Xc,ik(Tik )  =  C*  J  ri(t)  ■  s*k(t  -  Tik)  ej27rfcTlk  dt  .  (10) 

The  only  unknown  parameter  which  has  to  be  estimated  as  well  is  the  complex  radar  cross 
section  which  introduces  an  unknown  phase  shift.  The  integration  has  to  be  taken  over  the 
pulse  duration  of  the  transmitted  signal  sk{t). 

For  a  non-coherent  network,  where  no  phase  information  is  available,  Eq.  (10)  simplifies  to: 

Xnc,ik{rik)  =  a*lk  J  n(t)  ■  s*k(t  -rlk)dt  (11) 

The  unknown  complex  amplitudes  cqfe  have  no  impact  on  determining  the  time  delays  as 
they  are  independent  from  each  other. 

4.1  Position  dilution  of  precision  (PDOP) 

A  good  tool  for  visualizing  locating  accuracy  for  a  given  multistatic  radar  network  layout  is 
the  PDOP  mapping,  which  was  originated  with  launching  the  Loran-C  navigation  system  and 
came  into  much  wider  usage  with  GPS  [21]-[24].  The  dilution  of  precision  can  be  interpreted 
as  an  expression  which  describes  the  impact  of  the  positions  of  the  transmit  and  receive  nodes 
of  a  sensor  network  on  the  relationship  between  the  estimated  time  delay  and  the  localization 
errors.  Hence,  plots  of  PDOP  give  a  deep  insight  in  the  achievable  localization  accuracy  for 
a  given  set  multisite  radar  network  configuration. 

For  each  measurement  error  a  corresponding  dilution  of  precision  can  be  defined.  For  the 
two  dimensional  case,  where  X  =  (x,y),  it  is: 

J  +  a2 

PDOP  =  - y-  (12) 

or 

/  2 

TDOP  =  (13) 

or 

GDOP  =  VPDOP 2  +  TDOP 2  (14) 

where  <j2r  is  the  standard  deviation  of  the  range  estimation,  defined  by  standard  deviation 
of  the  time  delays  c2of ,  a2  and  a2  are  the  variances  of  the  localization  on  the  x  and  y 
axis,  respectively.  Furthermore  we  define  two  more  DOP’s,  one  for  the  horizontal  x-axis 
HxDOP  =  \fo2Joji  and  one  for  the  horizontal  y-axis  HyDOP  =  \foJJon. 


RTO-EN-SET-1 57(2010) 


6  - 11 


Aspects  of  Sensor  Networks 


ORGANIZATION 


In  Fig.  5  contour  maps  of  the  PDOP  for  a  non-coherent  (left  plot)  and  coherent  (plot  on  the 
right  side)  radar  network  are  shown.  The  multistatic  radar  network  consists  of  4  transmitters 
and  receivers,  equally  placed  on  a  circle.  Each  transmitter  emits  orthogonal  waveforms  and 
the  receivers  estimate  the  time  delay.  From  these  estimated  time  delays  the  target  position 
is  determined  by  solving  Eq.  (10)  or  for  the  non-coherent  case  Eq.  (11).  The  green  and  red 


Figure  5:  PDOP  contour  plot  with  4  transmitters  and  4  receivers  radar  network  symmet¬ 
rical  placed.  Left  figure  shows  the  accuracy  of  non-coherent  and  right  figure  the  coherent 
processing  case,  respectively.  (©IEEE  2009,  from  [20]) 


dots  in  the  PDOP  maps  represent  the  transmitters  and  receivers  of  the  radar  network.  For 
the  coherent  case  a  target  located  inside  the  network  circle  shows  lower  PDOP  values  than 
targets  located  outside.  In  comparison  to  the  non-coherent  case  the  accuracy  does  not  change 
much  for  target  inside  the  observation  circle  and  the  lowest  PDOP  value  is  obtained  at  the 
center.  Outside  the  radar  network  the  localization  errors  rises  dramatically  for  the  coherent 
case.  If  the  network  operates  non-coherently  the  best  PDOP  is  between  each  transmit-receive 
pair  and  degrades  smoothly  in  both  directions  from  the  baseline. 

In  the  second  example  the  transmitters  and  receivers  form  nearly  a  semicircle,  as  depicted 
in  Fig.  6  by  red  and  green  dots.  Clearly  visible  is  the  degradation  of  the  accuracy  of  mea¬ 
surements  for  the  coherent  case,  in  comparison  to  the  symmetrical  configuration.  These 
examples  show  that  high  localisation  accuracy  can  be  obtained  if  the  target  is  encircled  by 
the  nodes  of  the  radar  network.  Hence,  the  location  of  transmitters  and  receivers  has  to  be 
chosen  carefully  to  obtain  the  demanded  accuracy.  In  contrast  to  the  coherent  case  the  non¬ 
coherent  processing  does  not  show  significant  degradation  in  comparison  to  a  symmetrical 
setup. 

Maps  of  PDOP  represent  a  very  useful  tool  for  showing  the  accuracy  which  can  be  ob- 
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Figure  6:  PDOP  contour  plot  with  4  transmitters  and  4  receivers  asymmetrical  placed  radar 
network.  Left  figure  shows  the  non-coherent  case  and  on  the  right  hand  the  coherent  case. 
(©IEEE  2009,  from  [20]) 


tained  for  a  given  radar  constellation  or  for  choosing  the  best  node  locations  to  cover  a  given 
surveillance  area. 


5  Conclusion 

It  has  been  shown  that  multistatic  radar  network  performance  depends  not  only  on  radar  pa¬ 
rameters,  but  also  on  the  chosen  geometry. 

Overall  multistatic  radar  networks  have  many  advantages  and  can  increase  the  performance 
of  a  single  radar  unit  easily.  However,  these  benefits  come  at  the  cost  of  an  increased  com¬ 
plexity  and  a  high  demand  on  synchronising  all  nodes  in  the  network. 
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